THE DISTRIBUTION OF COSMOLOGICAL 7- RAY BURSTS 

Ehud Cohen and Tsvi Piran 
Racah Institute for Physics, The Hebrew University, Jerusalem, Israel 91904 

Received: 1994 Sept 5th; Accepted: 1995 Jan 25th 

Abstract 

We compare the burst distribution of the new (2B) BATSE catalogue to a cosmo- 
logical distribution. The observed distribution agrees well with a cosmological one, 
however, it is insensitive to cosmological parameters such as fi and A. The bursts are 
not necessarily standard candles and their luminosity can vary by up to a factor of 
ten. The maximal red shift, z max , of bursts longer than 2 sec is 2.1+ 1 7 (assuming no 
evolution). The present data is insufficient to determine maximal red shift, z max , of 
bursts shorter than 2 sec. 

Subject headings: gamma-rays: bursts - Cosmology: Observations 
1 Introduction 

The recent observations of the BATSE experiment on the COMPTON GAMMA RAY 
OBSERVATORY have suggested that 7-ray bursts (GRBs) originate from cosmological 
sources (Meegan, 1992) (see however, e.g. Podsiadlowski, Rees & Ruderman, 1994 for a 
different point of view). The count distribution, iV(C), in the BATSE IB catalogue agrees 
well with a cosmological distribution (Piran, 1992; Mao & Paczyhski, 1992; Dermer, 1992; 
Schmidt, 1992; Wickramasinghe et al. 1993, Mao, Narayan & Piran, 1994, denoted MNP). 
Additionally there is evidence (Norris, et al. 1994,1995, Nemiroff, et al. 1994), for the pre- 
dicted (Piran, 1992, Paczyhski, 1992) correlations between the duration, the strength and 
the hardness of the bursts (see however also Mitrofanov et al. 1994 and Fenimore 1994). 
With the release of the new (2B) BATSE catalogue (Meegan, et al. 1994) we analyze here 
the count distribution in terms of a cosmological model. 

We calculate the likelihood function to obtain the observed bursts' counts distribution, 
iV(C), from a variety of theoretical cosmological models. We show, first, that the observed 
count distribution in the 2B catalogue agrees with a cosmological distribution. The observed 
distribution is, however, insufficient to distinguish between different cosmological models. We 
demonstrate that the bursts are not necessarily standard candles and that their luminosity 
can vary by up to a factor of ten. We conclude by comparing our analysis with the results 
of Norris et al. (1994,1995) on the detection of cosmological time dilation and with the 
analysis of Fenimore et al. (1993) of the combined IB data and PVO data. 



1 



2 Basic Model 



We review here, briefly, the basic features of a cosmological source count distribution 
(see e.g. Weinberg, (1973) for a general population of sources and e.g. Piran, (1992) for 
GRBs). We consider a detector with a fixed energy range, AE } that operates for a fixed 
time span, At, and sources with a count spectrum: N(E) oc E~ a . The count rate at the 
detector from a source at a redshift z is: 

c(L z) - ^ l + z ) 2 - a m 

where L depends on the luminosity of the source, L, in the relevant energy range, AE } on 
the average energy E and on the observation time, At: L = L(AE)At/E. The luminosity 
distance is di(z). The number, iV(> C), of events with an observed count rate larger than 
C is: 

I' CO CZ(C,L) rp> f] r ( <y\ 

N(>C) = ^J Q n(L)dLj Q ^^ n {z)^jp-dz (2) 

where r p (z) is the proper distance to a redshift z and z(C, L) is obtained by inverting Eq. 1. 
The luminosity function is n(L) and the intrinsic rate (the number of events per unit proper 
volume and unit proper time at redshift, z) is n(L). The distribution, iV(C), depends on 
the cosmological parameters: and A, and on the source parameters: a, n(z) and n(L). 

We assume that the sources have a spectral index a = 1.5 (Schaefer et al. ; 1992). The 
resulting distribution is sensitive to the value of a and therefore we quote the parameters 
obtained for a = 2. The intrinsic evolution, n(z) and the luminosity function n(L) are 
unknown. We have chosen n(z) of the form: 

n(z) = n (l + z)-^ (3) 

Note that Piran (1992) uses n(t) = not 13 ' and for small z f3' —3/3/2. 
A standard candle distribution, 

n(L) = n 6(L - L ) } (4) 

is characterized by one parameter, z maX} the maximal z at the source is detected (obtained 
by inverting Eq. 1 for C = C m i n and L = Lq). The parameter n is determined from 
the analysis. We consider more general distributions that depend on two parameters. We 
examine a decreasing power law luminosity function with a minimal luminosity L : 

n(L) = [n (s-l)/L ](L/L )-% for L > L (5) 
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and an increasing power law with a maximal luminosity. These distributions are character- 
ized by z max ((L)) and by the power law index s. We also examine a maximally non - smooth 
distribution, a two standard candles distribution: 

n(L) = n 1 8(L-L 1 ) + n 2 8(L-L 2 ), (6) 

where the ratio rii/n 2 is chosen so that the number of bursts observed from each group is 
equal. This distribution is characterized by z max ((L)) and by the luminosity ratio, Li/L 2 . 

3 The Data 

The 2B catalog available in the public domain contains a list of all gamma-ray bursts 
that triggered the BATSE detectors up to 1993 March 9. The catalog contains the ratio 
Cmax/Cmin for 412 bursts, where C max is the maximum count rate and C m ,- n is the detection 
threshold and a duration table for 434 bursts (Meegan et al. 1994). BATSE's background 
varies as a function of time, so that the detection threshold C m ,- n does not remain constant. 
To obtain a regular sample we select a constant threshold C cut (C c „t( 1024ms) = 287 and 
Cc«t(64ms) = 72), and prune the data to contain only bursts which satisfy C m ,- n < C cut < 
Cmax (Petrosian, 1993; MNP). Following MNP we ignore the 256ms channel which seems to 
be more prone to irregularities. 

Kouveliotou et al. (1993) have shown that GRBs are bimodal in duration with short 
(^90 < 2 s) and long (6t 90 > 2 s) bursts. Since BATSE is more sensitive to the long bursts 
(MNP) it is meaningless to combine both sub-classes to one distribution and we have to 
separate the BATSE population to two sub-groups. In order not to exclude bursts which do 
not have a duration record we have used the relation (MNP): 

{C m a X /C cut ) W24: f > 1, for long bursts, 
(C max /C cut ) 64 X < 1, for short bursts, 
to add 15 bursts to our set. After excluding bursts which do not satisfy with all the previous 
requirements we have 238 long bursts, and 50 short bursts. 

In principle we would like to compare both the short and the long sub-populations to 
theoretical cosmological distributions. However, some of the short bursts are shorter than 
64ms. The count rate in this channel corresponds for these bursts to the total fluence of 
the burst and it underestimates the peak intensity of these events. This will distort the 
count distribution N(C). We exclude, therefore, the short bursts from this analysis and 
compare here only the long sub-class with the theoretical models. As far as we know all 
previous studies of the BATSE data ignored this problem and included the short bursts in 
their analysis. This study is the only one to date which proceeds to accomplish an accurate 
cosmological analysis and interpretation for LONG bursts only, using the 2B catalog. 
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4 Results 



The theoretical models depend on five unknown parameters; 0, A, z maX} s and f3 as well 
as on a. We calculate the likelihood function (see e.g. Press et al. , 1993) over this five 
dimensional parameter space and find the range of acceptable models (those whose likelihood 
function is not less than 1% of the maximal likelihood). We then proceed to perform a KS 
(Kolmogorov-Smirnov) test (see e.g. Press et al. 1993) to check whether the model is 
acceptable. We present in the following several two dimensional cuts through the parameter 
space. 

We begin with = 1 and A = and standard candles with no source evolution. For 
this simple model the likelihood function peaks at z max = 2.1 and the allowed range at a 
1% confidence level is: 1.4 < z max < 3.1 (for a spectral index a = 1.5). For a spectral index 
a = 2 the 1% confidence range is: zf^^ = l.bli'^. Figure 1 depicts the observed number 
count distribution for the long bursts (using the 1024ms channel) and the theoretical counts 
distributions for the allowed range of z max . 

Figure 2 depicts the contour lines of the likelihood function for the long bursts distri- 
bution in the (£l,z max ) plane. The likelihood function is practically independent of in the 
range: 0.1 < < 1. The likelihood function is also insensitive (figure not shown) to the 
cosmological constant A (in the range < A < 0.9, in units of the critical density). Following 
these results we are left only with the intrinsic parameters of the bursts distribution and we 
present in the rest of the paper results for cosmologies with = 1 and A = 0. 

Figure 3 describes the likelihood function in the (z maX} f3) plane for sources with varying 
intrinsic rate. The banana shaped contour lines show that a distribution with no intrinsic 
z dependence (/3 = 0) is equivalent to one with an increasing number of bursts with cosmo- 
logical time (/3 > 0) with a lower z max . This tendency saturates at high intrinsic evolution 
(large /3), for which the limiting z max does not go below .5 and at very high z maX} for 
which the limiting f3 does not decrease below -1.5. 

We examined three luminosity functions: a decreasing power-law with a minimal lumi- 
nosity, Eq. 5, an increasing power-law with a maximal luminosity and a double standard 
candles distribution, Eq. 6. Figure 4 depicts the likelihood function in the (s, z max ((L)) 
plane for luminosity functions of the form Eq. 5. We find that the likelihood function is 
quite insensitive to s as long as s ^ 2.2 (note that (L) diverges as s — > 2). The average 
luminosity with s = 2.2 is 5L 0} hence the luminosity varies, effectively by a factor of ten in 
this distribution. The likelihood function doesn't go below 70% of the maximal value for 
the whole range of power indices, s, for an increasing power law with an upper cutoff. This 
distribution is dominated by the strongest sources and the addition of weaker sources hardly 
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effects N(C). When we considered two standard candles sources, Eq. 6, the 10% likelihood 
is around a luminosity ratio of 10 between the sources. The KS test gives a probability of 
80% for a luminosity ratio of 14. These results demonstrate that the BATSE data alone do 
not require a narrow luminosity function. 

PVO has been detecting GRBs for about 14 years - but with a lower sensitivity than 
BATSE. The PVO count distribution, N(C') oc C -3 / 2 , corresponds to a homogeneous dis- 
tribution indicating that PVO is detecting only nearby GRBs whose distribution is not 
influenced by cosmological effects. In principle one should combine the BATSE and the 
PVO data. However, the PVO data are not available publicly at this stage. We can still 
incorporate the main PVO result in our analysis using synthetic data imitating the PVO 
data: a data set with the same number of bursts as observed by PVO, with an Eucleadian 
count distribution, N(C') oc C -3 / 2 , and a sensitivity limit larger by a factor of 20 than 
BATSE (Fenimore et al. (1993). We perform a maximum likelihood analysis (see Figure 
3) for the combined data and we find that the combined distribution h 0.87 (for 

standard candles with no source evolution). The allowed z max range at 1% confidence is 
0.65 < z max (PVO + BATSE) < 1.2. This reproduces the main conclusion of Fenimore et 
al. (1993) who combined the full (long and short) IB data with the PVO data. It is not 
known what is the fraction of long vs short PVO bursts. If we assume that the ratio of short 
to long bursts in PVO is the same as in BATSE, these results are slightly modified and we 
obtain: 0.55 < S max (PVO + BATSE) < 1.05. 

When we add luminosity functions (of the same kind described earlier) we find that the 
combined data require a narrow luminosity function. For example when we consider the two 
delta function distribution, the allowed luminosity ratio (at 1%) confidence is 2.5 compared 
to a factor of 10 in the BATSE data alone. Source evolution has the same effect as seen 
in the BATSE data - z max decreases when we consider sources whose intrinsic burst rate 
increases with cosmological time (see figure 3). 

The result that z max = 0.87 is the maximal red shift of the BATSE data should be 
compatible with the BATSE data on its own. However, a KS test rules out z max = 0.87 for 
the BATSE 2B data: The probability for z max = 0.87 of 2B data alone is less than 2 • 10~ 6 
(for a = 1.5). The 1% upper limit of PVO+BATSE (z max (PVO + BATSE) = 1.2) is ruled 
out at a 10~ 3 level. The discrepancy between the two data sets improves slightly when we 
include intrinsic source evolution (see Figure 3) but even with f3 = 2.5 the overlap between 
the likelihood functions of the two data sets is at a significance level of 0.3%. A possible 
alternative resolution is that this discrepancy arises from systematic errors in matching the 
sensitivity and the triggering algorithms of the two detectors. 
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5 Conclusions 

For standard candles with no source evolution the maximal red shift up to which the 
long bursts are observed is z max = 2.lt.o'^- With an estimated BATSE detection efficiency 
of 0.3, this corresponds to 2.3^°^ • 10~ 6 events per galaxy per year (for a galaxy density 
of 10~ 2 h 3 Mpc -3 ; Kirshner et. al, 1983). The rate per galaxy is independent of H and is 
only weakly dependent on 0. For = 1 and A = the typical energy of a burst whose 
observed fluence is Fr (in units of 10~ 7 ergs/cm 2 ) is 7^ • 10 5O -F7 ergs. These numbers vary 
slightly if the bursts have a wide luminosity function. The distance to the sources decreases 
and correspondingly the rate increases and the energy decreases if the spectral index is 2 
and not 1.5. The rate increases and the luminosity decreases if there is a positive evolution 
of the rate of bursts with cosmological time. 

Several luminosity functions with effective width of up to factor of ten fit the data. This 
is slightly wider than the allowed width (approximately a factor of six) obtained by Horack 
and Emslie (1994), who use a different statistical approach ("moment analysis") and whose 
data base is the BATSE IB catalogue. The two results are consistent in view of the differ- 
ent techniques, different data bases and different models for the luminosity functions used. 
This width is comparable to, and even wider than, the width of some observed luminosity 
distributions such as the luminosity functions of different types of supernovae. 

In agreement with Fenimore et al. (1993) we find that the combined BATSE data and 
the simulated PVO data corresponds to a cosmological distribution with z max 0.87. The 
combined data also require a narrow luminosity function. This result is, however, inconsistent 
with the 2B data on its own. A luminosity function or source evolution do not change this 
conclusion. This inconsistency suggests that there are systematic errors in matching the 
PVO and BATSE data. It might be rectified if either the sensitivity of PVO relative to 
BATSE has been overestimated, or if the PVO detection efficiency relative to BATSE is 
lower than expected. 

Norris et al. (1994,1995) found that the dimmest bursts are longer by a factor of 2.3 
compared to the bright ones (note however, that Mitrofanov et al. (1994) find no evidence 
for time dilation). With our canonical value of z max = 2.1 the bright bursts originate at 
^bright ~ 0.2. The corresponding expected ratio due to cosmological time dilation, 2.6, is in 
agreement with this measurement. If true, this results will support the cosmological model 
and imply that (i) intrinsic source evolution is insignificant and (ii) the low z max obtained 
form the combined PVO and BATSE data is rules out. Fenimore & Bloom (1994) find, 
on the other hand, that when other effects are included, this time dilation corresponds to 
Zmax > 6. In this case (ii) is still valid. However, (i) is invalid since z max > 6 requires a 
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significant negative intrinsic evolution: f3 ~ — 1.5 ± 0.3. 

The GRB count distribution is compatible with a cosmological one. Unlike preliminary 
hopes, (Trimble, 1994) we find that it is impossible to distinguish between different cosmo- 
logical models using the available GRB data. Firstly, N(C) is rather insensitive to and A. 
Additionally, like in other cosmological sources cosmological evolution can easily mask some 
cosmological effects. However, not all hope is lost. If GRBs result from neutron star merg- 
ers (Eichler et al. 1989), than detectors like LIGO and VIRGO could detect gravitational 
radiation signals from those events in coincidence with GRBs (Piran, 1992; Kochaneck & 
Piran, 1993). This would enable us to measure the distance to those objects (Schutz, 1986). 
This additional information, which is not available to any other cosmological objects, might 
enable us to determine cosmological parameters using the GRB distribution after all. 

We thank Tsafrir Kolatt, Jay Norris and Ramesh Narayan for helpful discussions and 
comments. This research was supported by a BRF grant to the Hebrew University and by 
a NASA grant NAG5-1904. 
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Figure 1: The normalized observed number counts distribution of the long bursts (with the 
1024msec channel) and the theoretical cosmological distribution. The theoretical distribution 
is shown for three values of z max : 2.1 (dotted line) 1.4 (long dashed line) and 3.1 (short dashed 
line) which correspond to the maximal likelihood value and the two 1% bounds. Also shown 
is a theoretical curve obtained for z max = 0.87 (dash-dot line: the best fit for the joint PVO 
and BATSE data). The descrepancy between this curve and the BATSE data is apparant. 
The curves are calculated for = 1, A = 0, a = — 1.5 and for standard candles with source 
evolution. 



Figure 2: Contour lines of likelihood function (levels 33%, 10%, 3.3% 1% etc.. of the max- 
imum, marked by +) in the (£l,z max ) plane for standard candles (Eq. 4), with no source 
evolution (/3 = 0) and a = 1.5. The cosmological constant A = 0. 



Figure 3: Gray scale map of the likelihood function (levels 33%, 10%, 3.3% 1% etc.. of the 
maximum) in the (f3 } z max ) plane for standard candles sources (Eq. 4) with a = 1.5 and 
"standard cosmology" = 1, A = for intrinsic evolution given by Eq. 3. Super imposed 
on these map is the similar likelihood function for the combined BATSE and PVO data 



Figure 4: Contour lines of likelihood function (levels 33%, 10%, 3.3% 1% etc.. of the maxi- 
mum, marked by +) in the (s, z max ((L)) plane for "standard cosmology" = 1, A = and 
source with no intrinsic evolution (/3 = 0) and a = 1.5 for a luminosity function given by 
Eq. 5. 
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